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SUMMARY 


The method of accounting for the effect of thickness on the lift 
distribution given in NACA Report 833 ^ 19^5 j has been modified. Compar- 
ison of lift distributions calculated by the methods of the present 
report and NACA Report 833 together with experimental distributions are 
presented. The results show that the modification improves the predic- 
tions of lift distributions for both angle of attack and camber. 


INTRODUCTION 


There are a number of methods available for the prediction of the 
lift distribution on given airfoil sections in incompressible flow 
(e.g., references 1 through 6 ) Of the methods cited, reference 3 alone 
presents a method which is relatively rapid and is applicable to 
arbitrary airfoils. The accuracy of the method of reference 3 , however, 
decreases with increasing amounts of section camber and/or thickness, 
particularly in the region of the airfoil nose. The purpose of the 
present report is to describe a modification to that part of the theory 
of reference 3 which deals with the lift distribution associated with 
angle of attack and/or camber. The modification is directed at improv- 
ing the accuracy of the method of reference 3 without decreasing its 
simplicity. No change is made in that part of the theory of reference 3 
dealing with the base profile. The discussion that follows presupposes 
a familiarity with the contents and method of reference 3 . 


METHOD 


The following relations for an airfoil of inf initesimal thic kne ss 
are presented in reference 3 (the notation used herein, unless otherwise 
noted, is the same as given in reference 3 and. is presented in the 
appendix) : 
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1. The velocity v induced at any point x Q on the camber line, 
due to all the vortices distributed along the camber line, is 


i P -a — ax 

V (*o) = 2n J x-x 0 (equation (l) of reference 3) (l) 

2. The bo undar y condition that the direction of the flow close to 
the camber line must be parallel to the surface of the camber line 
yields the relation 

= ^2- - a (equation (2) of reference 3) (2) 

V 0 dx 

3. The distribution of vorticity is assumed to be 

= 2 7 0 (aq 1 cot I + E A n sin n0^ 

(from equation (4) of reference 3) (3) 

4. The lift at a point along the chord is 

cLqT 

Q p = (equation (l4) of reference 3) (4) 

5. The total lift is 

° L = So p7 ° “dx* dx ^ 

The coefficient A 0 ‘ in equation (3) is shown to vary only with 
angle of attack, while the coefficients A n are independent of the 
angle of attack and are functions of the mean camber-line shape only. 
Thus' the part of the lift distribution which varies with angle of attack, 
called the additional— lift distribution, and the part of the lift dis- 
tribution which is a function of the mean camber— line shape only, called 
the basic— lift distribution, can be treated separately. 

As stated in reference 3 > the preceding equations were derived from 
a theory for thin airfoils where the local velocity at each elemental 
vortex along the camber line was taken to be free— stream velocity V 0 . 
Reference 3 suggests that, for the airfoil of finite thickness, a more 
accurate prediction of the lift distribution would, be obtained if the 
local velocity at each vortex station were taken to be the local velo- 
city Vf at the surface of the base profile at the same chordwise 
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station. Reference 3 approximates the effect of the local velocity by 
adjusting the local lift at each point by the ratio of local to free— 
stream velocity. That is. 



(equation (28) of reference 3) (6) 


Equations ( 4 ) and (5) then become, for the airfoil of finite thickness, 

' PV f 




P = 


dx 


( 7 ) 


and 



( 8 ) 


The approximation given by equation (6) led to improved accuracy of the 
lift-distribution calculation. No attempt was made in reference 3 to 
improve further the accuracy by .solving for a new vorticity distribution 
using the new boundary condition which would result if were substi- 

tuted for V 0 in equation (2). 


The need for further improvement in the prediction of pressure 
distributions, with particular emphasis on accounting for the effects of 
thickness — or in effect to account for the fact that Vf becomes 
markedly different from V 0 — has made a re-examination of the problem 
necessary. In an attempt to calculate more accurately the effects of 
thickness on the chordwise load distribution, the boundary condition 
given by equation (2) is altered so that it becomes 


v _ <tyc _ 
V f dx 


( 9 ) 


With this new boundary condition a solution for the vorticity distribu- 
tions can now be made in a manner similar to that used in reference 3* 


The expression for the vorticity distribution is again assumed to 
be given by equation (3), but now equation (9) is used to define the 
coefficients of the series. From equations ( 1 ), (3), and ( 9 ) it can be 
shown that 

— (— 2 - - a) = -rA 0 * + 2 An cos n 9 (10) 

V Q V <1* / 1 

The coefficients are given by 



( 11 ) 
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and 

‘ f I Q % (iGT 008 ”» 119 < 12 > 


The coefficients can each he separated into two parts, one part a func- 
tion only of .the camber— line shape, and the other part a function only 
of the angle of attack. Following previous work, the vorticity distri- 
bution, ascribed to the camber-line shape, is that existing when there is 
zero vorticity at the nose of the airfoil and is called the basic dis- 
tribution. The angle of attack at which this distribution exists is 
called the ideal angle < 4 . Therefore, the total angle a is written 
as the sum of the ideal and additional angles of attack, as follows 


a = a i + a*. 

Equations (11) and (12) can be written as 




and 


- °i) cos 1,0 39 - 1 1 j- 


cos n0 d0 (14) 


It i$ now possible to write the expressions for the coefficients corre- 
sponding to basic and additional vorticity. For basic vorticity these 
are 


and 


N> = I / - “0 008 1,9 39 


(15) 


(16) 


The coefficients for the additional vorticity are 



and 



cos n 9 d 0 


(17) 

( 18 ) 


Thus the basic— and additional— vorticity distributions for the airfoil 
finite thickness can be expressed as 
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and 



(19) 

( 20 ) 


Basic— Lift 'Distribution 


•Since the basic vorticity is zero at the nose of the airfoil it. ia 
required that the coefficient Aq^' be equal to zero. This requirement 

fixes the value of aq which is necessary for the evaluation of the 
coefficients A n ^. Setting equation (15) for An * equal to zero, 
that is. 



the value of aq 


is obtained as 


a i 



Tfo dx 


de 


d 9 
Vo ' 


( 21 ) 


The basic— vorticity distribution is now written as 



.= 2V Q L A-,-, sin nQ 
u 1 n b 


( 22 ) 


The basic lift is derived from equations ( 7 ), (16), and (22) as 

P *o = 4 (vo} I I Q f “ a i ) do3 ”0 sin nflo d 0 ( 23 ) 


which can also be written 




V f dy ct 
Vo ** 


cot 



(24) 
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where 



and. 

d ^c h d ^c 
dx dx 1 

In a similar manner, the camber required for a specified "basic lift 
is- obtained as 



where 



The accuracy of the basic— lift distribution calculated by use of 
equation (24) will be shown to be an improvement over the distribution 
calculated by the expression given' in reference 3* 


Additional— Lift Distribution 
♦ 

The additional lift can be derived from "equations (7), (17), (18), 
and (20). At the point 9 0 the additional lift is 


= 4 


o 




Vf 

V 0 


cos tid sin n0 o 



(26) 


Use of equation (26) to express the additional— lift distribution results 
in a loss in computational simplicity as compared with the expression 
for the additional lift as given in reference 3. In addition, a compari- 
son with experimental results indicates a loss in accuracy. However, a 
simple and adequate method for calculating the additional— lift distri- 
bution has been found by using a different approach. Throughout the 
remainder of the report it will be assumed that the basic lift is 
evaluated by the method just presented, and given by equation (24), and 
the additional— lift distribution is evaluated by the following, method. 
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The additional— vorticity distribution of the infinitesimally thin 
airfoil, given by 

= 2V 0 A 0 ' cot | (27) 

a 

is adjusted in proportion to the change in local velocity from V Q to 
Vf at each vortex station. The resulting expression for the additional— 
vorticity distribution of the finite— thickness airfoil is 



= 2V f A 0 * cot | 


(28) 


The additional— lift distribution is obtained by direct substitution 
of this expression for the vorticity distribution in the expressions for 
the local additional lift and the total additional— lift coefficient. 

The latter two expressions are, for the finite— thickness airfoil, respec- 
tively. 


( 

P V-p \ rl Y J 
1 o \ ai ^a ( 


(29) 


and 





qc 



( 30 ) 


The expression that is obtained for the additional lift per unit total 
additional— lift coefficient is 




cot 

2 


+ cos 9 ) d0 


( 31 ) 


It may be noted that equation (31) can also be written in a form similar 
to equation (56) of reference 3, that is. 





( 32 ) 
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rather than 



(equation (55) of reference 3) 


The accuracy of additional— lift distributions calculated by use' of 
equation (31) will be indicated later by comparing them with experimental 
lift distributions. 


METHOD OF APPLICATION 


In applying the method of this report for the calculation of the 
lift distribution on an airfoil section, it is necessary to know the 
velocity distribution over the base profile. There are several ways to 
obtain the required base— profile velocity distribution. Reference 4 
and a number of subsequent MCA reports, such as reference 7, contain 
velocity distributions for a large number of MCA base profiles. A 
convenient method of calculating the velocity distribution on base 'pro- 
files is presented in reference 3* After the base— profile velocity 
distribution has been determined, the basic— and additional— lift 
distribution can be calculated by use of equations (24) and (31), respec- 
tively. 


The total lift at a point for a given value of t ota 1- s ec t i on— li ft 
coefficient Cj is 



where 



(equation (57) of reference 3) 


r 


P. d(x/c) 


(33) 


Equation (24), for the basic lift, can be solved numerically by either 
the method given in the appendix of reference 3 or by the method of 
reference 8. If, as is generally the case, the analytical relation 


between — and 9 is not known, 

T 

v o 

lift, can be solved graphically. 


equation (31) ^ for the additional 


The inverse problem of calculating an airfoil section which will 
have a specified lift distribution can be solved in the manner outlined 
in reference 3- Equation (25), given herein, is used in determining the 
camber— line slope and can be solved by either the numerical method given 
in the appendix of reference 3 or by the method of reference 8. 
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EVALUATION OF METHOD 


An evaluation of the additional—, basic— , and total— lift distribu- 
tions calculated by the method of this report is provided by comparing 
them with lift distributions obtained by the method of reference 3 and 
by experiment. Comparisons of additional— lift, basic— lift, and total- 
lift distributions are shown in figures 1, 2, and 3> respectively. The 
experimental data in figure 1(b) were obtained from reference 9 and the 
experimental data in the remainder of the figures were obtained from 
unpublished data from the Ames 7— by 10-foot wind tunnel. For the com- 
parisons of basic-lift distributions (fig. 2),' experimental data were 
not available for the ideal angle of attack as given by either the 
method of this report or by reference 3* ■ Therefore, the available 
• experimental data were interpolated to the ideal angles of attack as 
given by equation (21) and by reference 3 • 

From the. comparisons shown in the figures, it is evident that the 
lift distributions calculated by equations (2k) and (31) of this report 
are in better agreement with the experimental results shown than are the 
distributions calculated by the method of reference 3; it will be noted 
in partic ular that the agreement is better in the region near the air- 
foil leading edge. Further, although the method of this report is based 
on the assumption of small camber, as in the case of thin airfoil theory, 
the calculated basic— lift distributions for the sections with relatively 
large camber (fig. 2) are in good agreement with experiment. 


CONCLUDING REMARKS 


Computation of lift distributions by the method of this report 
gives results which are in closer agreement with experimental distribu- 
tions for the data examined than those obtained by use of the airfoil 
theory given in reference 3. The method of this report, in effect, 
modifies the method of reference 3 by adjusting the vorticity, given by 
thin-airfoil theory, for ‘the local velocity induced by the airfoil- 
thickness distribution. 


Ames Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Dec. 12, 1950. 
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APPENDIX 

NOTATION 

coefficients of the trigonometric series for distribution of 
vorticity 

section chord, feet 
section lift coefficient 
total lift on the section, pounds 
local lift coefficient at any point x 
free— stream dynamic pressure,, pounds per square foot 
induced velocity at any point x, feet per second 
free— stream velocity, feet per second 

local velocity at the surface of the base profile, feet per 
second 

distance along chord, feet 
mean line ordinate, feet 
angle of attack of the chord line, degrees 
ideal angle of attack of the chord line, degrees 

rate of change of circulation, feet per second 
trigonometric coordinate , cos' 
density of free-stream air, slugs per cubic foot 

Prescript 

indicates application to airfoil of zero thickness 




o 



Subscripts 


indicates application to additional loading 

indicates application to basic loading 

indicates application to specific point along chord 
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Calculated by method of this report 

Calculated by method of reference 3 

O Experiment 


O 7 .2 .3 .4 .5 .6 .7 .8 .9 


Chord station, £ 




(a) NACA 63,012 

Figure l - Comparison of the calculated and experimental 
additional -lift distributions. 



Additional lift per unit total additional-lift coefficient. 
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(b) NACA 0018 
Figure /. - Concluded. 



Basic lift, P h Basic lift 
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(a) NACA 64A8/0, a =0.8 mod. 



(b) NACA 0010-63, Cj. - 0.8, a =0.8 mod. 


Figure 2. - Comparison of calculated and experimental 
basic - lift distributions . 





Total lift, P Total lift 
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(a) NACA 64A8I0, a =0.8 mod., c t = 0.93. 



(b) NACA 00/0-63, C,. =0.8, a =0.8 mod; c t =0.65 

Figure 3.- Comparison of calculated and experimental 
total- lift distributions. 


Total lift, 
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Chord station, -jf 


(c) NACA 00/0-63, c t . -0.8, a - 0.6 mod.; c t - 0.90 
Figure 3. - Concluded. 
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